Two-proton radioactivity and three-body decay. V. Improved momentum 

distributions. 
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Nowadays quantum-mechanical theory allows one to reliably calculate the processes of 2p ra- 
dioactivity (true three-body decays) and the corresponding energy and angular correlations up to 
distances of the order of 10"^ fm. However, the precision of modern experiments has now become suf- 
ficient to indicate some deficiency of the predicted theoretical distributions. In this paper we discuss 
the extrapolation along the classical trajectories as a method to improve the convergence of the the- 
oretical energy and angular correlations at very large distances (of the order of atomic distances), 
where only the long-range Coulomb forces are still operating. The precision of this approach is 
demonstrated using the "exactly" solvable semianalytical models with simplified three-body Hamil- 
tonians. It is also demonstrated that for heavy 2p emitters, the 2p decay momentum distributions 
can be sensitive to the effect of the screening by atomic electrons. We compare theoretical results 
with available experimental data. 

PACS numbers: 21.45.-v, 21.60.Gx, 23.50.-l-z 



I. INTRODUCTION 

Two-proton radioactivity is the most recently discov- 
ered radioactive decay mode of nuclei and it is a very 
actively developing field. There were 42 years between 
the prediction [ij and discovery [H, Q of 2p radioactivity 
and, subsequentially, seven years latter we have several 
well studied examples. A number of experiments per- 
formed in the last 2-3 years can be characterized as key 
for the field. In particular, correlations in the 2)d decays 
have been measured recently in ^Be [4], ^^Ne [5j, ^^Mg 
0, Hi, ''^Fe 0, and ^'^Ag [1| providing qualitatively new 
information about the 2p decays. With correlation in- 
formation becoming available, the 2p decay studies are 
now turning into a field of research where precise infor- 
mation about structure and continuum dynamics can be 
obtained. It is clear that our ability to extract useful in- 
formation from correlations is directly dependent on how 
well we understand the propagation of particles in the 
long-range three-body Coulomb field. 

From a theoretical point of view, true two-proton decay 
(2p radioactivity) is an exclusively quantum-mechanical 
phenomenon, which has no analogue in classical physics. 
It is expected to be widely spread along the proton drip 
line with Z < 50 due to peculiarities of the pairing in- 
teraction. A consistent quantum-mechanical theory of 
two-proton radioactivity and "democratic" three-body 
decays of the coulombic nuclear systems has been de- 
veloped in the series of papers j9l4l3| , which we continue 
here, and has been applied to different physical cases in 
Refs. [1, [T3 - [l7j . The complete momentum correlations 
for the decay of a non-aligned three-body system can 



be described by two parameters. These parameters are 
chosen in this and our previous studies as the energy dis- 
tribution parameter e between any two of the particles 
and the angle 9k between the Jacobi momenta: 

e^E^/Et , cos(0fc) = (k, • ^)/(A:, fcj,) , (1) 
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where Ai are mass numbers of the constituents, M is a 
nucleon mass, and Et = Q2p is a two-proton decay en- 
ergy. For two-proton emitters these parameters can be 
constructed in two "irreducible" Jacobi systems, called 
"T" and "Y" , see Fig. [H The detailed definition of the 
Jacobi coordinates can be found in Ref . \4!\ ■ The complete 
correlation pictures for two- pro ton decay were, for the 
first time, calculated in Ref. Various aspects of the 
correlations between the decay products have been dis- 
cussed in the theoretical works of Refs. [3, 0, [HI [H, 

In ^Be and ^^Fe, the complete correlation pictures for 
2p decay were recently obtained experimentally 0, 0]- 
Moreover, the precision of these experimental results is 
now sufficient to show a deficiency in certain aspects 
of the predicted momentum distributions in the case of 
heavy 2p emitters [13] . It was already understood in Ref. 
[lOl] that this deficiency is connected to the limited radial 
range of the calculations and the approximate nature of 
the boundary conditions employed for the treatment of 
the three-body Coulomb asymptotic. 
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FIG. 1. Independent "T" and "Y" Jacobi systems for the 
core+A'^+TV three-body system in coordinate and momentum 
spaces. There are "planar" cases where both the coordinates 
and the momenta belong to the same plane. 



The classical extrapolation (CE) of momentum distri- 
butions was suggested in Ref. [ll| as a simple way to es- 
timate the possible influence of the "residual" Coulomb 
interaction. The basic idea is that, at small distances, 
particles are propagated by quantum-mechanical equa- 
tions providing the three-body wave function (WF) ^'3^' 
with outgoing asymptotic. At some sufficiently large dis- 
tance, the WFs are converted into "events" with definite 
coordinates and momenta by a Monte-Carlo (MC) pro- 
cedure. However, at that time (in 2002 the 2p decay of 
''^Fe was just discovered with statistics of the order of 
ten events [1, Q) the need to improve this aspect of our 
calculations was assigned to the remote future and no 
detailed studies were performed. Now it seems that the 
development of the field has achieved the stage where the 
need to improve this aspect of our approach has become 
evident. 

In this work we discuss the method of classical ex- 
trapolation in detail, demonstrate its reliability by ap- 
plication to exactly-solvable three-body models with a 
simplified Hamiltonian, and consider three "key" cases 
(^Be, ^^Mg, and ^^Fe) covering a broad range of possible 
charges, masses, and structures for the 2p emitters. 

Natural system of units with ?i = c = 1 is used in this 
work. 



II. APPROXIMATE BOUNDARY CONDITIONS 

In this section we sketch the methods used to construct 
the approximate boundary conditions 10] and outline 
existing problems. The asymptotic form of the three- 
body potentials in the hyperspherical harmonics (HH) 
method is 



Vr]K'y,K''y' 



where multiindex {Kj} = {K, S^lx-,ly, s^} is a com- 
plete set of quantum numbers. The matrix Uk-i,K'-i' 
arises due to contributions from the short-range nuclear 
forces, and NK-f,K'-t' > are some integer numbers. The 
effective contribution of the short-range forces decreases 
as p^^ or faster in hypersherical space. The diagonal 
centrifugal term depends on the "effective angular mo- 
mentum" £ = K + 3/2. Coulomb pairwise potentials 
generate the long-range part of the hyperspherical po- 
tentials behaving as p~^. From the technical side, the 
three-body Coulomb interaction causes problems due to 
long-range channel coupling (nonzero nondiagonal "Som- 
merfeld parameters" 'qK-i,K''y') that does not allow one 
to decouple the HH equations on the asymptotic. To 
deal with this problem, the finite-size potential matrix 
(in truncated hyperspherical basis) can be diagonalized 
with respect to the long-range term by the orthogonal 
transform V = A^VA: 
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(3) 

This potential includes nondiagonal "centrifugal" terms 
CK-y,K'-y' and, to achieve the asymptotic in the diagonal- 
ized representation, we still need to go very far in p value, 
where the terms ~ become negligible compared to 
those with ~ . At such p values, the hyperradial part 
of the asymptotic solution with pure outgoing nature can 
be constructed in the form 



(+) 



ip) ^ ^K',,K'Y {Gco iVK'Y , p) + iFc„ {llK'Y ,p)), 



K'j' 



^3^^ = p-'^' E x^Kiip)JK,in,) . (4) 

The functions F and G are the ordinary regular and 
irregular Coulomb functions. Hyperspherical harmon- 
ics J'K-y are functions of the 5-dimensional "solid angle" 
^5 — {SpT^xi^y}- Here and f2j, are ordinary solid 
angles of the Jacobi vectors X and Y [see Eq. [5] and 
tan(6'p) = y/Mx/MyX/Y. The value £0 should be larger 
than 3/2, but otherwise does not seem to be particularly 
important. The WFs x^"*"^ provide the necessary bound- 
ary conditions for the decay problem. 

The proposed boundary conditions are exact on the 
truncated hyperspherical basis at a hypersphere of very 
large radius. However on a practical level, these two re- 
quirements contradict each other: the movement further 
in radius requires the increase of the basis size; a larger 
basis size may require a larger radius. Therefore, at some 
point, the further radial propagation of the solution (with 
fixed basis size) leads to a deterioration of its quality. For 
■*^Fe with the decay energy of 1.154 MeV and the basis 
size of iCrnax = 20, radii between 500 and 2000 fm are 
needed to get reasonable solutions. 

There exists an analytical asymptotic of the three- 
body Coulomb problem (a so called "Redmond- 
Merkuriev" asymptotic [l^ 13), which is presumably 
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applicable to the true three-body decay. Practical appli- 
cation of this asymptotic is technically complicated and 
it seems that there exists a very limited experience in us- 
ing such an asymptotic. At the moment we are going to 
avoid these complexities and to demonstrate that there 
exists a simple and practical way to treat the problem. 



III. EXTRAPOLATION ALONG CLASSICAL 
TRAJECTORIES 

To perform a classical extrapolation of the quantum- 
mechanical result, we need to switch from a WF to clas- 
sical trajectories. This should be made at some closed 
surface around the decay region. The procedure becomes 
especially simple if the whole surface is located in the re- 
gion of classically allowed motion. Then the flux vectors 
at the surface can provide initial conditions for classical 
trajectories. 

When using the HH coordinates there is only one 
variable p, which has a dimension of length [the 6- 
dimensional flux can be calculated for different p values, 
see Eq. ([S])]. Therefore, it is natural in this approach to 
select the hypersphere with a large radius Pmax as such a 
surface. We will see later that tiny regions on the hyper- 
sphere with a large radius, where the pairwise distances 
appear to be small, do not lead to problems as the WFs 
in these regions are strongly suppressed. This happens 
due to the energy conditions defining the true 2p decay: 
there are no long-living states in either pair of the three 
constituents and the strong Coulomb repulsion rapidly 
"expels" particles from the regions where they are close 
to each other. 

A less evident, but important requirement is that the 
hyperradius /Omax is large enough that the typical dis- 
tances between each pair of particles significantly exceeds 
the typical quantum coherence length (the "corpuscular" 
aspect of the problem is then far prevailing over the pos- 
sible wave effects). This is a complicated issue and in 
each case an acceptable minimal value of Pmax should be 
defined by numerical experiment. 

The classical trajectories formed at this hypersphere 
Pmax are propagated to distances poxt 3> Pmax at which 
the momentum distributions are stabilized (what this ex- 
actly means we will see below). After this, the momen- 
tum distributions are reconstructed from the set of tra- 
jectories. 

The pairwise distances, the Jacobi vectors, and the 
hyperradius are connected by the following relations 



ri2 = X , r23 = Y - ciX , rgi = Y + C2X , 



„2 _ MA2 ^ (Ai + ^2)^3 



Ai + A2 A1 + A2 + A3 
ci = Ai/{Ai + A2) , C2 = A2/{Ai -f A2) . (5) 

In the definition of the hyperradius p, particle A3 should 
be a heavy core if X and Y are defined in the "T" Jacobi 
system and either Ai or A2 should be a core in the "Y" 



Jacobi system (see also Fig.[l]for the numbering conven- 
tion) . 

The Newton equations of the motion for the Jacobi 
vectors are used to avoid the extra degrees of the freedom 
connected to the center-of-mass motion: 
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The particular choice of the form of Eqs. © ("T" or "Y" 
Jacobi system) or the numerical precision in solving this 
system are not practical obstacles for getting the correct 
classical trajectories. 

The initial conditions for these equations are defined 
on the hypersphere of the maximal radius achieved in the 
quantum-mechanical calculations: 

{Pmax,f^^'-'}^{X(0),Y(0)}, 

{Upn...M;^)MPn...Mp)} ^ {x(o),Y(o)}, (7) 

where f2p is a randomly generated 5-dimensional hyper- 
angle selected by the MC procedure according to the WF 
density l^'s'''^ P at p = Pmax- The fiux associated with the 
Jacobi vectors is defined in an ordinary way: 



-Im 
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(+) 



(8) 



In the quantum- mechanical model of the three-body 
decays |9l-[l3j. the total flux j through the hypersphere 
p = Pmax define the width 



(9) 



where A'' is normalization of the WF ^'3~'~'* in the in- 
ternal region. The momentum distribution (distribu- 
tion density) is found as the differential of the fiux 
dj /[de dcos{9k)], see Eq. In this work we com- 

pare the quantum-mechanical distributions calculated at 
P = Pmax (called below "without classical extrapolation" 
or "initial") with distributions obtained by classical ex- 
trapolation to p = pext ( "with classical extrapolation" or 
"final"). 



A. Treatment of spins 

It is implied above that the flux is averaged over the 
initial spin states and summed over the final spin states. 
Therefore the components of the WF ^'3''''' with different 
total spin S can be considered as different "particles" 
whose contributions to the total momentum distribution 
should be added incoherently. 

In general, three particles (or two Jacobi vectors) de- 
fine a plane. Within this plane, the set of 6 equations ^ 
can be reduced to 4 equations. However, the momentum 
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vectors do not necessarily belong to this plane. It is evi- 
dent that the geometry of the problem remains planar in 
the case of zero angular momenta of the X and Y sub- 
systems (this situation is shown in Fig. [IJ. For nonzero 
angular momenta, some additional considerations are re- 
quired. 

Let us consider the flux field induced by the ordinary 
two-body WF with I ^ 0. For to = 0, the flux is purely 
radial as the angular part of the WF Yim is real (flux is 
an imaginary part of the gradient matrix element). For 
purely radial flux, the classical angular momentum asso- 
ciated with the particular trajectory is zero (radius and 
momentum vectors are coUinear). This can be seen as 
a source of a confusion as the quantum-mechanical mo- 
mentum of the WF and the classical momentum of the 
selected trajectory are explicitly different. The answer 
seems to be that the classical characteristic of the trajec- 
tory should be related to average corresponding charac- 
teristic of the WF. 

In the three-body case, the ground-state WFs typi- 
cally have two major components: the dominating L = 
component and an "admixture" L = 1 component. We 
imply here that a spin-zero core is considered; the two 
spin 1/2 protons can then be coupled into the total spins 
5 or 5 = 1. The L = component of the WF is 
formed by terms with angular momenta in the subsystem 
Ix = ly It is easy to check that the angular part of this 
WF [Y;^ (8)11^] 00 is real and thus the classical angular mo- 
mentum associated with any trajectory induced by this 
WF is zero. The decay in this case is planar (we mean 
that for any generated event, a plane can be selected in 
which the coordinate vectors and the momentum vectors 
of all three particles are simultaneously located). 

It is more complicated when the L = 1 component is 
considered. It is possible to demonstrate that for the 
[Yi^ (g) Yi^jiM component of the WF with M = 0, the 
configuration of the classical momenta is planar, while 
for M ^ 0, the planes formed by the three radii and by 
the three momenta does not coincide. However according 
to the Wigner-Eckart theorem, in order to define the ob- 
servables it is sufficient to calculate the matrix elements 
for only one projection and the rest are reconstructed by 
the angular momentum algebra. Therefore, it seems suf- 
ficient to calculate the distributions for M = (planar 
case calculations are especially simple), while the distri- 
butions for M = ±1 should be the same. 



IV. TEST CASES OF SOLVABLE 
SEMIANALYTICAL MODELS 

In Ref. a semianalytical model was developed 

which allows one to treat exactly the asymptotic be- 
haviour of the three-body Coulomb WF for certain sim- 
plified three-body Hamiltonians. The basic idea of the 
model is that instead of the real three-body Hamiltonian 

H3^Tx + Ty + Vuiru) + T^23(r23) + ^^3i(r3i) , (10) 



we use the model Hamiltonian depending not on pairwise 
vectors but on the Jacobi vectors X and Y 

H3^Tx + Ty + Vx(K) + Vy{Y)+V3{p), (11) 

The three-body potential V3{p) is used in this work has 
the Woods-Saxon form 

Vsip) = V,' (1 + exp [{p - po)/ap])"' , (12) 

PO = V2 1.2(Acorc + l)'^\ (13) 

with a small value of the diffuseness parameter ap = 
0.4 fm. The depth of this potential is used to control 
the decay energy of the system. The potentials Vx and 
Vy contains the nuclear and the Coulomb contributions. 
The Coulomb potential of the homogeneously-charged 
sphere with a radius Vgph is used. The nuclear parts are 
described by Woods-Saxon formfactors with radii taken 
from systematics. 

In conjunction with this simplified Hamiltonian of Eq. 
pil) . we can introduce an auxiliary Hamiltonian 

H3 = Tx + Ty + Vx(Xi2) + VyiY23), (14) 

for which the Green's function can be constructed in an- 
alytical form 

(XY, X'Y') = —J dEx (X, X') Gi,+j (Y, Y') , 

(15) 

where Et is the total decay energy. Ex — sEt, and Ey = 
(1 — £)Et are the energies of the Jacobi subsystems. The 
above two-body Green's functions correspond to the X 
and Y subhamiltonians of H3. Based on Eq. (|15p . the 
width and the energy distribution for the system defined 
by the Hamiltonian of Eq. (fTTj) can be obtained from 

dV dj 8 MxMy 2 f,^. 

where dj/de is a differential of the flux at the asymptotic. 
For a particular set of quantum numbers Ix, ly, the ampli- 
tudes ^(e) are defined via the scattering eigenfunctions 
ipi. of subhamiltonians of (|14l) : 

/>oo />oo 

A(e)=/ dX dYipiJkx{e)X)^i^{ky{e)Y) 

xV3{p)ipLU,s{X,Y). (17) 

The WF (p]:^i^i^g{X,Y) is the quasistationary eigenfunc- 
tion of pT|) . deduced in a three-body hyperspherical ap- 
proach. The particular choice of the boundary conditions 
for this WF (for sufficiently large radius of the "box" ) is 
not important in the model. The quasistationary WF is 
normalized to unity in the internal region, which gives 
the identity dT/de = dj/de in Eq. (fT6)) . 

The results obtained in this model are quoted be- 
low as "exact" as they do not suffer from any conver- 
gence/stability issues. In Sections IIV Al and IIVBI we 
will use models with different simplified Hamiltonians to 
test the classical-extrapolation procedure in the case of 
the ^^Mg-ground-state (g.s.) decay and only after that 
we will turn to more realistic situations. 
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FIG. 2. Classical trajectories for ^ Mg in a direct decay 
model (pmax = 1000 fm, Et = 0.75 MeV). Panel (b) shows 
one selected trajectory on a large scale. The dotted line in 
panel (b) corresponds to the scale of the panel (a). 
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A. Direct-decay model 

The Hamiltonian of Eq. (jlip constructed in the "Y" 
Jacobi system corresponds to some physically well justi- 
fied approximations. Namely, (i) we neglect the proton- 
proton interaction and (ii) for one of the core-proton po- 
tentials we use the Jacobi Y variable instead of the rel- 
ative distance. The later assumption becomes correct in 
the limit of an infinitely-heavy core and thus should work 
well for heavy 2p emitters. 

Let us consider the "Y" system, where the subsystem 
{core-|-proton} is taken as an effective particle lying on 
the X coordinate as shown Fig. [T] 

In this case, we include both pairwise interactions V^f"*^' 
and Vy^'^^. The system with such a composition of poten- 
tials in the "Y" system was labeled as "Two final-state 
interactions" in Ref. [l2j . 

For the i^^Mg g.s., we assume the pure d-wave structure 
= = 2 in this model. The nuclear Woods-Saxon 
potential was used with the radius 

ro = 1.2(^ + 1)'/' (19) 

and diffuseness a = 0.65 fm. The depth of the potentials 
was adjusted to give an energy of 1.3 MeV for the ground- 
state resonance in -i^Na Q and the Coulomb potential of 
the charged sphere with radius 

r,ph = (1.2^1/3)2 + 0.82 (20) 

was used. In the above expressions, one should substi- 
tuted A — A2va the X subsystem and A — A2 + l'v[\ the 
Y subsystem. In this model, we obtained the half-life of 
Ti/2 = 58 ps (corresponding to F = 7.9 x 10"" MeV) 
which is in qualitative agreement with the experimental 
value for ^^Mg (T1/2 = 4 ps [6]). 

The radial convergence of the energy distribution e in 
this model for some classical trajectories is illustrated in 



FIG. 3. (Color online) Energy distribution for ^^Mg with 
different pmax values without (upper panel) and with (lower 
panel) classical extrapolation in the direct decay model. The 
calculations are performed with Et = 0.75 MeV, /Cmax ~ 20, 
and lower panel with pext = 40000 fm. Gray curve shows the 
exact result of Eq. (|16p (the same for both panels). 

Fig. [21 The trend of the CE is to make the energy distri- 
bution narrower. The visual stability of the distributions 
is achieved at distances of about poxt ~ 7000 fm [Fig. 
[5fa)]. On a larger scale, a certain drift of the trajectories 
can be seen up to much larger distances [Fig. [2jb)]. 

The effect of the CE on the energy distributions is 
demonstrated in Fig. |31 The energy distributions have 
a characteristic bell shape. The upper panel shows 
the energy distributions calculated with the quantum- 
mechanical three-body model fl3| for different pmax val- 
ues. The calculated result tends towards the "exact" 
result of Eq. (|16p . shown by the gray curves. However, 
this convergence is very slow and some discrepancy re- 
mains even for the largest available Pmax- The lower 
panel shows the distributions obtained with the classi- 
cal extrapolation. These distributions are clearly wrong 
for Pmax ^ 500 fm. However for larger pmax, they stabi- 
lize and reproduce the results of the solvable model [Eq. 
P^ ] within the width of the curve. 



B. "Diproton" model 

The word "diproton" in parenthesis is the name of this 
model as it is different from the diproton model typi- 
cally used in the literature. The diproton correlation 
in our model is not introduced statically (which means 
"by hand") but is treated dynamically. In Ref. [llj, we 
have demonstrated that when introduced appropriately 
for the configurations with lowest possible angular mo- 
menta in the subsystems, the diproton model can provide 
only a very small value for the 2p width. For decays of 
the higher-Z configurations, like [p^] or [d?] for 0"^ states, 
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FIG. 4. Classical trajectories for ^''Mg in the "diproton" 
model (pniax = 1000 fm, Et = 0.75 MeV). Panel (b) shows 
one selected trajectory on a large scale. The dotted line in 
panel (b) corresponds to the scale of panel (a). 

this model overestimates the width. Therefore, it is not 
apphcable in practice, in contrast to widespread beliefs. 

In this work we apply the diproton model, not for re- 
alistic estimates, but for testing purposes. The diproton 
model gives very sharp energy distributions focused at 
small p-p energies. So, we use it to determine whether 
the CE procedure works for conditions of strong kine- 
matical focusing. 

In the diproton model, Eqs. ([TTt - pT| are used in the 
"T" system, where the core {^3,^3} interacts with the 
two protons as if they were an effective particle {Ai + 
A2,Zi + Z2}. The Coulomb potential of the simplified 
Hamiltonian can be written in the form 

ycoui ^ aZiZ2 ajZi + ^2)^001.0 , . 

X Y ■ ^ ' 

Note, that this is a model with only one nuclear pair- 
wise interaction V^^^'^{X) in the p-p channel (the second 
interaction can be put to zero) and therefore the model 
is called "One final-state interaction" in Ref. [l^l- The 
proton-proton nuclear potential for an s-wave is taken as 
a single Gaussian 

y(r) =-K,exp[-(r/ro)2], (22) 

with Vb = — 31 MeV and ro = 1.8 fm reproducing the 
low-energy s — nucleon-nucleon phase shifts. The 
Coulomb potential of the charged sphere with radius 

rsph - {l.2Aliiy + ^{1.2x2^/3)^ (23) 

is used in the Y coordinate. The half-life of ^^Mg ob- 
tained in this model is T1/2 = 0.39 ps (corresponding to 
r = 1.2 X 10-9 MeV). 

The radial convergence of the energy e in this model 
for some classical trajectories is illustrated in Fig. H) The 
trend of the CE is for the trajectories to drift towards the 
more narrow "diproton" peak in the energy spectrum. 
The convergence trend is analogous to the direct-decay 
model with several thousand fm required for a reason- 
able stabilization and more than a hundred thousand fm 
required for complete stability. 




FIG. 5. (Color online) Energy distributions in ^^Mg for 
the "diproton" model without (upper panel) and with (lower 
panel) classical extrapolation. Calculation results are shown 
for different pmax values. The calculations are performed 
with Et = 0.75 MeV, A'max = 14 and for lower panel with 

Pcxt — 

10^ fm. Gray curves show the "exact" result of Eq. 
(|16|) (the same for both panels). 

The effect of the CE on the energy distribution 
is demonstrated in Fig. [5j The case appears to 
be completely analogous to the direct-decay model. 
The upper panel shows the energy distributions calcu- 
lated within our three-body hyperspherical quantum- 
mechanical approach for different Pmax values. The 
quantum-mechanical results tend towards the "exact" re- 
sult (|16p . but only very slowly. The distributions pro- 
vided by the classical extrapolation (see the lower panel 
in Fig. [S]) contain artifacts for pmax ^ 500 fm, but for 
larger pmax, they stabilize and reproduce the result of the 
solvable model Eq. (|16l) within the width of the curve. 



C. Brief conclusions 

Before we continue studies of realistic cases, let us out- 
line what we can conclude on the basis of the exactly- 
solvable models with simplified Hamiltonians. 

(i) The quantum-mechanical calculations performed 
for pmax of a few thousand fm give energy distribu- 
tions which have visible deviations from the "exact" 
results obtained in the semianalytical model. The 
extrapolated distributions practically coincide with 
the "exact" ones. 

(ii) The CE provides decent results only if the start- 
ing point for the extrapolation is sufficiently large. 
Pragmatically, this means that the classical trajec- 
tories in the kinematical space {e,cos{9k)} should 
be quite short. The same should be true in the 
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conjugated coordinate space. It can be expected 
that the criterion of a successful transition from 
quantum to classical calculation is that the classi- 
cal variation of a position in some space should be 
smaller than corresponding coherence length. 

(iii) Distances of tens of thousands of fm's are needed 
to achieve complete stabilization of classical trajec- 
tories in practice. Some very minor drift continues 
after that, reflecting the long-range nature of the 
Coulomb interaction. However, it is evident that 
distances of 10^ fm are already atomic scale dis- 
tances and the nuclear Coulomb effects should be 
suppressed for larger distances due to some form of 
electron screening. 

Near perfect convergence of the extrapolated distri- 
butions to those calculated in the exact semianalytical 
models with the simplified Hamiltonians is not a proof 
that the procedure should work perfectly in the case of 
a complete three-body Hamiltonian. However, it is very 
encouraging and we can expect that the quality of conver- 
gence in the realistic case will be very similar, as the kine- 
matical conditions for the decay in the simplified models 
are chosen to be the same as in the realistic cases. 



REALISTIC THREE-BODY CASES 




FIG. 6. Contour maps of the distribution density on the kine- 
matical plane {e, cos{dk)} for ''Be in "T" (upper row) and "Y" 
(lower row) Jacobi coordinate systems without (left panels, 
"init.") and with (right panels, "fin.") classical extrapola- 
tion. 



For the models with simplified Hamiltonians, we 
demonstrated only the energy distributions (angular dis- 
tributions are trivial) and only in one Jacobi system 
(the one in which the particular semianalytical model 
is formulated). Conversion of the distribution into the 
other Jacobi system in this case does not provide an ad- 
ditional information. For realistic calculations we will 
demonstrate complete correlation pictures (on kinemat- 
ical {e,cos{9k)} plane) simultaneously in both "T" and 
"Y" Jacobi systems. It should be understood that cor- 
relation pictures in "T" and "Y" Jacobi systems are 
just different representations of the same physical phe- 
nomenon. Conversion between these distributions is triv- 
ial. Nevertheless, we systematically demonstrate both of 
them simultaneously as each representation allows one to 
reveal different aspects of the correlations (see Ref. 
as example). 

A. Decay of the "^Be 

Very precise complete correlation data have recently 
been obtained for ^Be in Ref. The detailed theoretical 
studies of 2p decay of the ^Be 0+ ground state have been 
carried out in this work and compared to the experimen- 
tal data. The dynamical range of around Pmax ~ 1000 
fm used in these calculations was estimated in Ref. 
as sufficient for essentially complete convergence of the 
momentum distributions. A very nice agreement between 
theory and experiment was found in this work. We would 



like to check here whether the conclusions obtained in 
Ref. Q can be influenced by a more careful treatment of 
the momentum distributions. 

The classical trajectories for ^Be in the kinematical 
space are all very short. Only the trajectories corre- 
sponding to small initial inter-particle distances [e ^ 0.5, 
cos{9k) ^ ±1 in "T" system] have noticeable lengths. 
The complete correlation densities without and with ex- 
trapolation are shown in Fig. |5] (this is the calculation 
with potential set P2 from Ref. [J], which was found to 
be the optimal choice in that work). The distributions 
are very similar except for the aforementioned regions 
of small initial inter-particle distances. A closer look at 
these regions is provided in the inclusive distributions in 
Fig- El The maximal effect can be found at small e values 
(corresponding to the lowest relative-energy motion be- 
tween two of the particles) or for the angular distribution 
in the middle energy bins around cos{0k) ~ ±1 (in the 
"T" system) and cos{9k) ~ — 1 (in the "Y" system). 

Comparisons with experimental angular distributions 
[3] are shown in Fig. [5] The theoretical curves here are 
visibly distorted (relative to Fig. [T]) as the comparison 
is based on the full MC simulation of the experimental 
setup [3] , which takes into account the effects of the ex- 
perimental bias and resolution. The effect of the classical 
extrapolation is at the limit of the experimental sensitiv- 
ity. Quantitatively the x'^/'^ values without extrapola- 
tion are 1.17 (in the "T" system) and 1.14 (in the "Y" 
system). The same values with extrapolation are found 
as 1.20 and 1.16, respectively. This is a little worse, but 
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FIG. 7. (Color online) Inclusive energy and angular distribu- 
tions for ®Be in "T" (upper row) and "Y" (lower row) Jacobi 
coordinate systems without (solid curves) and with (dashed 
curves) classical extrapolation. Black lines show the total 
distribution and the color coded lines show the inclusive dis- 
tributions for certain energy and angular bins (described in 
the legends). 



not really significant. On the other hand, there seems to 
be a minor improvement of the agreement for the parts of 
the middle energy bins mentioned in the previous para- 
graph. 

The properties of the ^Be continuum are actively in- 
vestigated now. The new higher precision experiments 
have been performed recently at NSCL (Michigan State 
University) and at Flerov laboratory ( JINR, Dubna, Rus- 
sia). The expected precision of these experiments would 
make the improvement of the theoretical distributions 
introduced in this work a necessary part of the data in- 
terpretation. 



B. Decay of the ^''Mg 

A systematic view of the classical trajectories on the 
kinematical plane for ^^Mg is given in Fig. [3] The 
"lengths" of the trajectories here are significant: typi- 
cally around 10 — 15% of the kinematical variable range 
thus making the CE procedure necessary for quantitative 
calculations of the momentum distributions. 

An improvement of the momentum distribution due to 
classical extrapolation is demonstrated in Fig. [TU] for the 
complete momentum distributions and in Fig. [Tljfor the 
inclusive ones. It can be seen that the angular distri- 
butions in the "T" and the energy distributions in the 
"Y" Jacobi systems are the most sensitive to the extrap- 
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FIG. 8. (Color online) Comparison of experimental (data 
points 01) and predicted (curves) cos(0fc) distributions in the 
"T" (left) and "Y" (right) Jacobi systems for the indicated 
gates on £ = E^/ Et parameter. The solid (green) and dashed 
(red) curves correspond to the three-body calculations with- 
out and with classical extrapolation. The effect of the detector 
bias and resolution is included. 



olation. The effect of the extrapolation on the distri- 
butions in certain energy and angular bins can be very 
large. The energy distribution in the "T" system is only 
slightly modified by the CE, but it is interesting to note 
that for very small e values (where p-p Coulomb inter- 
action is expected to be most active) the extrapolated 
distribution is visibly suppressed. 

Unfortunately the available experimental data on the 
momentum distributions in ^^Mg do not provide com- 
plete distributions, but provide distributions projected 
on a plane (perpendicular to the incident beam axis). 
Such distributions integrated over one variable have lost 
some information and can be more complicated to inter- 
pret. 



C. Decay of the ^^Fe 

The ^^Fe nucleus is the heaviest 2p emitter studied so 
far and the effect of the CE is the largest, see Fig. [121 

Radial stabilization of the values e and cos(6'fe) for one 
selected trajectory is demonstrated in Fig. [13] (this trajec- 
tory is shown in the gray ellipse in Fig. [T^ . The trajec- 
tories are well "converged" by about (3 — 4) x 10* fm but 
some drift continues up to much larger p values. In real 
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FIG. 9. (Color online) Classical trajectories on the kinemat- 
ical plane {e,cos{6k)} for ^^Mg in the Jacobi "T" system. 
Starting points (larger black dots) correspond to pmax ~ 1000 
fm. Dots in the curves correspond to pcxt equal 1300, 2000, 
3500, and 10^ fm. The red dot at the axis cos{Ok) ~ corre- 
sponds to a stationary point, see the discussion in Sec. IVI Cl 
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FIG. 11. (Color online) Inclusive energy and angular distri- 
butions for ^^Mg in "T" (upper row) and "Y" (lower row) 
Jacobi coordinate systems without (solid curves) and with 
(dashed curves) classical extrapolation. Black lines show the 
total distribution and the color coded lines show the inclusive 
distributions for certain energy and angular bins (described 
in the legends). 
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FIG. 10. Contour maps of the distribution density on the 
kinematical plane {e,cos(Sfc)} for ^^Mg in "T" (upper row) 
and "Y" (lower row) Jacobi coordinate systems without (left 
panels, "init.") and with (right panels, "fin.") classical ex- 
trapolation. 



experimental situations, this slow drift can be suppressed 
by electron screening which is discussed separately in Sec. 

Em 

The improvement of the momentum distributions due 
to the classical extrapolation for ''^Fe is demonstrated in 
Fig. HH for the complete momentum distributions and in 
Fig. [15] for the inclusive ones. The most impressive mod- 
ifications are for the e distribution in the "Y" system and 



FIG. 12. (Color online) Classical trajectories on the kine- 
matical plane {e,cos(&fc)} for *^Fe in the Jacobi "T" system, 
Et ~ 1.154 MeV. Starting points (larger black dots) corre- 
spond to Pmax = 1000 fm. Dots in the curves correspond to 
pext equal 1400, 2200, 4000, and 10^ fm. The red dot at the 
ajcis cos(^fc) = corresponds to a stationary point, see the 
discussion in Sec. IVI CI 



for the cos(6'fe) distribution in the "T" system. As far as 
these distributions have bell shapes, centered at (or close 
to) the center of the kinematical range, we can charac- 
terize them in terms of the full width at half maximum 
(FWHM) . Classical extrapolation decreases this value by 
about 30% for cos{9k) in "T" system and by about 10% 
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FIG. 13. Radial stabilization of the values e [panel (a)] and 
cos(0fc) [panel (b)] with p in the case of screened Coulomb 
potential (solid curves) and in the case of nuclear Coulomb 
potential only (dashed curves) for one selected trajectory in 
^^Fe (see Fig. [HI)- Et = 1.154 MeV, p^ax = 1000 fm. 




FIG. 14. Contour maps of the distribution density on the 
kinematical plane {e,cos{9k)} for *^Fe in "T" (upper row) 
and "Y" (lower row) Jacobi coordinate systems without (left 
panels, "init.") and with (right panels, "fin.") classical ex- 
trapolation. 



for e in "Y" system. This effect is sufSciently large to be 
already observable at the current level of the experimen- 
tal precision. 

The experimental distribution for ^'^Fe |7] has quite 
low statistics (150 events) and therefore it is far from 
being smooth, see Fig. [16] (a). To make a visual compar- 
ison with theoretical calculations possible, we produce a 
"smooth" representation of this data based on the ex- 
perimental uncertainties. The raw experimental data 
measured in Q by an optical time projection chamber 
consists of the energies and the polar angles of the two 
protons and the azimuthal angle between the projections 
of the two protons' momenta on the cathode plane of the 
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FIG. 15. (Color online) Inclusive energy and angular distri- 
butions for *^Fe in "T" (left column) and "Y" (right column) 
Jacobi coordinate systems without (solid curves) and with 
(dashed curves) classical extrapolation. Black lines show the 
total distribution and the color coded lines show the inclusive 
distributions for certain energy and angular bins (described 
in the legends). 



chamber. Each parameter for each event has a value (and 
its uncertainty) defined individually by a complex itera- 
tive fitting procedure. Instead of each event, we generate 
an event distribution based on the stochastic Gaussian 
variation of each parameter within its uncertainty range. 
So instead of one point in the kinematic space we get 
a kind of a "probability cloud" . The result of this pro- 
cedure is shown in Fig. [TCj (b). This procedure is not 
a cure for small statistics, but for small statistics and 
large experimental uncertainties we think it is a prefer- 
able presentation as it incorporates information about 
the distortions caused by the measurement procedure in 
a consistent and visible way. 

The experimental data are compared with inclusive 
theoretical distributions sensitive to the classical extrap- 
olation in Fig. [T7| In this plot, theoretical results were 
treated by the procedure which is maximally close to the 
experimental treatment of the data: (i) for the "theoreti- 
cal event" the nearest experimental event in the space of 
parameters {Epi, Ep2,0i,O2, \4>2 — 0i|} was defined, (ii) 
spherical coordinates for protons from the "theoretical 
event" were distributed according to the errors of the 
nearest experimental event, (iii) the momentum of the 
core was reconstructed and the total energy of the "dis- 
torted" theoretical event was renormalized to correspond 
exactly to the experimental one, and (iv) a new location 
in the kinematical plane {£,cos(0fc)} was defined. The 
effect of the experimental resolution is a roughly 25% 
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FIG. 16. Experimental distribution density in Fe in "T" 
Jacobi system. Original distribution from [3] is shown in the 
panel (a) as a histogram. "Smooth" version of this distribu- 
tion taking into account experimental errors is shown in the 
panel (b) as a contour plot. 



VI. DISCUSSION 
A. Classical motion 

It is important to note that large hyperradii are used to 
start the classical extrapolation procedure. Specifically 
for true 2p decay with such large hyperradii, practically 
the whole WF resides is in the classically allowed region 
(probability to find the system in the classically forbidden 
region is very small). For example, for the ""^Fe calcula- 
tion with hyperradius pmax — 1000 fm and MC genera- 
tion of 10^ events it is typical that not a single event is 
generated which is situated in the classically forbidden 
region. This fact confirms the validity of the choice of a 
hypersphere as the surface at which the switching from 
quantum-mechanical to classical methods is performed. 
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FIG. 17. (Color online) Inclusive angular (a) and energy (b) 
distributions for *^Fe in "T" and "Y" Jacobi coordinate sys- 
tems without ("init.") and with ("fin.") classical extrapo- 
lation compared to the experimental data. The results with 
experimental resolution of Q and without it are explained in 
legends, which are the same for both panels. 



increase of FWHM for the e distribution and an 18% in- 
crease of FWHM for the cos(6'fc) distribution (see Fig. 
[T7| . It can also be seen in Fig. [T7] that the theoreti- 
cal results with classical extrapolation are in quantita- 
tive agreement with the experiment, while without CE 
they are not completely consistent with the data. So, we 
have appreciable experimental evidence that the long- 
range treatment of the momentum distributions (namely 
CE) is necessary for heavy 2p emitters. 



B. Electron screening 

The discussion of the *^Fe case can provide an illus- 
trative example here. So far, the decay process of ^^Fe 
with the half-life of 2.6 ms [7| was measured in gas (or 
solid state) detectors. This means that at the moment of 
decay, '*^Fe has completely recovered electron shell. The 
Bohr radius for '*^Fe is 



1 



ao 



m,eaZ 



= 2035 fm , 



(24) 



where Z = Zi is total charge of ^^Fe. Therefore, we 
can expect that the screening effect of the inner most 
electrons becomes observable at about 2000 fm. Clas- 
sical trajectories for ^^Fe in kinematical space are well 
stabilized by 10^ fm, but there is a minor drift up to 
much larger distances. It is clear that some effect of the 
electron screening on the momentum distributions can be 
expected. 

The binding energy of all electrons estimated as inde- 
pendent particles is '^^me{Z /2niaY [rii is a principal 
quantum number of the shell), which gives 52.3 and 47.8 
keV for ''^Fe and ^"^Cr, respectively. So, when ''^Fe emits 
two protons at least two electrons should be ejected car- 
rying away 4.5 keV of energy. The estimated velocities 
of protons with energies around 0.5 MeV and electrons 
with energies around 1 keV are 0.033 and 0.063. These 
velocities are comparable, which means that the 2p decay 
of atomic ^^Fe would be accompanied by a strong recon- 
struction of atomic structures having the same timescale. 
It is reasonable therefore to make estimates of a screening 
with the ^^Fe electron density, but only for 24 electrons. 
This will somehow account for the effect of the electron 
shell disintegration during the 2p decay of ^^Fe and pro- 
vide a nuclear plus atomic Coulomb potential tending to 
zero at infinity. 

The electron density used for the screening calculations 
and the potentials obtained are shown in Fig. [151 One 
can see already that at 2000 fm, the full (V^^i + Koui) 
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FIG. 18. Panel (a) shows electron density for 24 lowest elec- 
tron shells in "^^Fe (normalized for integration over dr). Left 
axis of panel (b) shows proton potential for Coulomb interac- 
tion of nucleus, electron shell, and their difference (screened 
potential). The ratio of the screened potential to the nuclear 
is shown by gray curve opposite the right axis. 



Coulomb potential is noticeably reduced due to the 
screening compared to nuclear Coulomb potential (the 
reduction factor is 0.8). At 7000 fm the reduction factor 
is 0.5 and it tends to zero at 30000 fm. 

The radial stabilization of the values e and cos(0fc) in 
the screening case compared to the purely nuclear case 
is shown in Fig. [13] for one selected trajectory. It can be 
seen that in the screening case, the trajectory stabilizes 
at p ~ (3 — 4) X lO'' fm. In the purely nuclear case, the 
minor drift of the trajectory continues to much larger p 
values. The calculations shovir that in the "T" system, 
the screening effect is largest for the variable cos(0fe). 
It is typically at the level of 0.6% of the absolute value 
of this variable and, for pmax — 1000 fm, it typically 
accounts for 3 — 4% of the CE effect. For an effect which 
is 0.6% at the absolute scale it is difficult to speculate 
about its observability just now: its scale is comparable 
to the widths of the lines in our plots. However, if we 
think about it as an effect of the atomic surrounding 
on nuclear decay properties, then such a value can be 
considered as an impressive one. 

It should be noted that the existence of the screen- 
ing effects is the subject of the experimental technique 
employed. For example, the 2p decay in ^^Mg was stud- 
ied in the decay- in- flight experiment Ref. In this 
experiment the ^^Mg g.s. was populated by the neutron 
knockout from the relativistic beam of the completely 
stripped ^°Mg ions. The resulting ^^Mg is also com- 
pletely stripped and can hardly pick up any electrons 
before the decay. Therefore, in spite of a long lifetime 
(^1/2 = 4 ps, which is much longer than typical recombi- 
nation time), screening in this experiment will have dif- 
ferent character, compared to the case discussed above 
for 45Fe. 



{e, cos(6'fc)} plane in the "T" system. For such points, the 
classical trajectories in this plane have zero length. For 
the degenerate situation e = 1, the stationary behavior is 
trivial; this situation is not very interesting as the phase 
space for such configurations in the quantum-mechanical 
problem tends to zero. However there exist nondegen- 
erate stationary points, and which for two-proton decay 
with a heavy core {^3, Z3} are found as: 



(1 + A3/2)(Z3M3)^/^ 
2A3 (Z3M3)2/3 



=) = o. 



(25) 



These stationary points are defined by the condition that 
the force acting on each particle is always directed exactly 
along the line connecting that particle with the center of 
mass of the whole three-body system. Such stationary 
points should exist for any two-body potential with the 
same power dependence on radius V(r) ~ r" for each 
pair of the particles. The values of e equal to 0.497, 0.382, 
0.379 are found for ^Be, ^^Mg, and ''^Fe, respectively, by 
Eq. (f25|) as well as by a direct calculations using Eq. ([6]) . 
It is clear that the solution, which is a stationary one 
in the {e,cos(0/c)} plane, is an analogue of the Lagrange 
solution in celestial mechanics (with reservation that we 
are dealing here with repulsive 1 jr potentials) . 

The multi-cluster decays of nuclear systems has been 
qualitatively studied in Ref. [l^l- In this work a quasi- 
classical approach was used, based on the classical self- 
similar solutions of the few-body Coulomb problem. The 
stationary point discussed above represents such a self- 
similar solution in our specific case. It was concluded 
in Ref. 20] that "Three-cluster configuration asymptot- 
ically approaches to an expanding self-similar triangle 
whose sides obey the (M/Z)^/'^ rule." This statement 
is probably not completely correct. It can be seen from 
Figs. l9l and [T2] that there is a trend for classical trajecto- 
ries to tend somehow towards the stationary point which 
corresponds to a self-similar solution. This trend leads to 
certain systematic modifications of the momentum distri- 
butions by the long-range Coulomb interaction. However 
as we have seen in this work, the whole picture is more 
complex. The total distributions occupy broad regions 
of the kinematical plane. They are determined mainly 
by the internal structure of the three-body system and 
the decay dynamics under the barrier, than by the long- 
range Coulomb interaction outside of the barrier. Clas- 
sical trajectories originating on the hypersphere of large 
radius each converge to its own final position which, in 
a general case, could have nothing in common with a 
stationary point. 



VII. CONCLUSION 



C. Self-similar solutions 

Prom Figs. 12] and jT^] it is possible to see that there 
exist so called "stationary points" in the kinematical 



In this work we discuss the extrapolation along the 
classical trajectories as a method to improve the momen- 
tum distributions for radioactive 2p decay (true three- 
body decay). The proposed method provides near per- 
fect description of the distributions in the test cases of 
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simplified three-body Hamiltonians. In the case of real 
three-body Coulomb interactions considerable quantita- 
tive effects on the distributions are observed. In the case 
of the lightest 2p-emitter ^Be this effect is minor, but in 
the heavier 2p-emitters (^^Mg and "^^Fe) the improvement 
is essential for the precise description of the distributions. 

It should be emphasized that some aspects of the mo- 
mentum distributions for 2p decays are sensitive to the 
long-range three-body Coulomb interaction, while the 
others are absolutely insensitive. Namely, the angular 
distribution in the Jacobi "T" system and the energy 
distribution in the Jacobi "Y" system are considerably 
modified by the classical extrapolation. Two other inclu- 
sive distributions (the energy distribution in the Jacobi 
"T" system and the angular distribution in the Jacobi 
"Y" system) are essentially not influenced by the classi- 
cal extrapolation. Therefore the long-range part of the 
three-body Coulomb does not practically change the in- 
formation about the internal structure of the decaying 
system which is contained in the latter distributions. 

Attention should be paid to the huge range which is re- 
quired both for the extrapolation range 10^ fm) and 
for the starting point of the classical procedure (~ 10"^ 
fm) at typical decay conditions. The classical procedure 
is applicable only for distances above 500 — 1500 fm (in 
p variable) for the considered set of 2p emitters (which 
is actually quite representative). The intermediate dis- 
tances from 30 — 100 fm (where the protons come from 



under the Coulomb barrier) to around 1000 fm should to 
be treated quantum mechanically to obtain decent results 
from the classical extrapolation. 

We have shown that the electron screening can have a 
sizable effect on the momentum distribution in the 2p de- 
cay of atomic ^^Fe. So, the 2p radioactivity belongs to a 
rare class of nuclear phenomena, which exist on the bor- 
derline with atomic phenomena. There exist examples of 
weak radioactive decay modification induced by atomic 
electrons (e.g. due to the energy conditions making /3~ 
decay possible only into bound electron states ^21] or due 
to the hyperfine effect [l^). We think that the sizeable 
sensitivity of the radioactive decay via particle emission 
due to a modification of the potential barrier properties 
in the atomic environment is demonstrated in our work 
for the first time. 
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